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ABSTRACT 

Context. Spitzer & wide-field surveys and followup capabilities have allowed a new breakthrough in mid-IR spectroscopy up to red- 
shifts > 2, especially for 24/im detected sources. 

Aims. We want to study the mid-infrared properties and the starburst and AGN contributions, of 24 fim sources at z ~ 2, through 
analysis of mid-infrared spectra combined with millimeter, radio, and infrared photometry. Mid-infrared spectroscopy allows us to 
recover accurate redshifts. 

Methods. A complete sample of 16 Spteer-selected sources (ULIRGs) believed to be starbursts at z ~ 2 ("5.8/um-peakers") was 
selected in the (0.5 deg 2 ) J1064+56 SWIRE Lockman Hole field ("Lockman-North"). These sources have SiAiim > 0.5 mJy, a stellar 
emission peak redshifted to 5.8 pm, and ry ega > 23. The entire sample was observed with the low resolution units of the Spitzer/IRS 
infrared spectrograph. These sources have 1.2 mm observations with IRAM 30m/MAMBO and very deep 20 cm observations from 
the VLA. Nine of our sources also benefit from 350 pm observation and detection from CSO/SHARC-II. All these data were jointly 
analyzed. 

Results. The entire sample shows good quality IRS spectra dominated by strong PAH features. The main PAH features at 6.2, 7.7, 8.6, 
and 11.3 //mhave high S/N average luminosities of 2.90±0.31, 10.38± 1.09, 3.62±0.27, and 2.29±0.26x 10 10 L G , respectively. Thanks 
to their PAH spectra, we derived accurate redshifts spanning from 1.750 to 2.284. The average of these redshifts is 2.017±0.038. This 
result confirms that the selection criteria of "5.8 /jm-peakers" associated with a strong detection at 24 pm are reliable to select sources 
at z ~ 2. We have analyzed the different correlations between PAH emission and infrared, millimeter, and radio emission. Practically 
all our sources are strongly dominated by starburst emission, with only one source showing an important AGN contribution. We have 
also defined two subsamples based on the equivalent width at 7.7 pm to investigate AGN contributions. 

Conclusions. Our sample contains strong starbursts and represents a particularly 24 /jm-bright class of SMGs. The very good cor- 
relation between PAH and far-IR luminosities is now confirmed in high-z starburst ULIRGs. These sources show a small AGN 
contribution to the mid-IR, around ~20% or less in most cases. 

Key words. Galaxies: high-redshift - Galaxies: starburst - Galaxies: active - Infrared: galaxies - Submillimeter - Techniques: 
spectroscopic 



1. Introduction 

The mid-infrared (MIR) regime is known to provide very 
rich diagnostics of the interstellar medium (ISM) of galaxies. 
Hot dust from starbursts and especially AGN contributes to con- 
tinuum emission. Dust also produces various absorption fea- 
tures, most spectacularly from silicates. However, the MIR spec- 
tra of star-forming galaxies are dominated by strong emission 



* Colour figures and the Appendix are only available in the electronic 
form via http://www.edpsciences.org 



features attributed to polycyclic aromatic hydrocarbons (PAHs). 
These large molecules, typically a few hundred carbon atoms 
in the ISM, are the carriers of the series of prominent emission 
bands observed between 3 and 19/zm, especially between 6 and 
12fim ( |Leger & Puget 1984||Aflamandola et al. 1985| l. PAHs ef- 
ficiently absorb UV and optical photons from young stars and 
reemit the energy mostly in the mid-infrared. 

After their discovery in the local ISM and nearby star- 
forming regions, PAH features have proved to be common, not 
only in various regions of the Milky Way, but in other galaxies 
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ranging from local systems to high-z starbursts (e.g., |Tielens| 
2008 ). PAH features were the focus of a good fraction of the 
MIR spectroscopy programs of the infrared space observatories 
ISO and Spitzer. Important key projects of both missions were 
devoted to MIR spectroscopy of large samples of galaxies. The 
ISO Key Project on Nearby Galaxies (Helou et al. 2000 1 showed 



an amazing similarity of PAH spectra in various galaxy envi- 
ronments. More detailed studies with the SINGS spectral map- 
ping program ( Kennicutt et al. 2003 1 using the Spitzer Infrared 



Spectrograph (IRS: Houck et al. 20041 provided further evi 



dence of interband PAH feature strength variations, especially in 



the presence of weak AGN (Smit h etal. 2007| >. 

The high sensitivity of SpitzerflRS allowed the extension of 
PAH spectroscopy to z ^ 2 for samples of hundreds of galax- 
ies that could be observed in the ~ 5-35 pm spectral range. This 
led to the measurement of a large number of redshifts of dust- 
enshrouded Spitzer sources at 1.5 < z < 2.5 in the so-called 
"redshift desert" of optical spectroscopy. These MIR spectra al- 
low one to distinguish power-law AGN-dominated spectra from 
PAH starburst-dominated spectra and composite systems. This 
long series of high-z IRS studies has thus established the cri- 
teria for the identification of various types of sources from 
Spitzer 3. 5-24 pm photometry ( |Fadda et al. 2010 Houck et al. 
2W5\ | Yan et al. 2005] |2007| JSajina et al. 2007| |Weedman 



et al. 2006[ |Dasyra et al. 2009| |Hernan-Caballero et al. 2009 
Bertincourt et al. 2009| l. Particular types of sources can be se- 



lected using various criteria: dust-enshrouded AGN from strong 
24 pm fluxes and/or power-law IRAC 3. 6-8. 0pm spectra (Dey 



et al. 2008| |Dey & The N dwfs/MIPS Collaboration 2009 ), star- 
bursts from the presence of a visible redshifted 1.6 pm stellar 
bump (|Farrah et al. 2006} |2008| [Lo nsdale et al. "20091 |Huang 



et al. 20091 jDesai et al. 20091), Type 1 or Type 2 QSOs ( |Lutz 



et al. 2008{ |Martinez-Sansigre et al. 2008 1, and submillimeter- 
selected galaxies (SMGs: Valiante et al. 2007| IMenendez- 



2009 Po pe et al. 2 008 Co ppin et al.] 



Delmes tre et al. 2007 
2010] !. 

Indeed, strong PAH emission is so characteristic of starbursts 
that the strength of PAH features may be used as a star forma- 
tion rate indicator. This is justified by the good correlation that 
has been found between the luminosity of PAH features (mostly 
at 6.2 and 7.7 pm) and the total IR or far-IR luminosity, which 
in the absence of an AGN is taken to be a good tracer of star- 
forming activity ( |Kennicutt 19 98 ). Such a correlation has been 
checked for starburst galaxies at both low and high redshift (e.g., 



Brandl et al. 2006; Pope et al. 2008 Menendez-Delmestre et al 
2009[ ). It is the origin of the correlation found between the ISO 



7.0pm- or Spitzer 8.0/Utn-band flux and the star formation rate 
(e.g., |Elbaz et al. 2002[ ~ 



Wu et al. 2006 Alonso-Herrero et al. 



2006 Farrah et al. 2007 1 in low-z luminous infrared galaxies. 
It is similarly related to the use of the single-band 24 pm lumi- 
nosity for determining star formation rates of galaxies at high 
redshift, especially at z ~ 2 (e.g., |Bavouzet et al. 2008} IRieke 
et al. 2009||Elbaz 2010| Fiolet et al. in prep.). 

Such rich diagnostics of the properties of the ISM and star 
formation at high redshift provided by mid-IR PAH spectral fea- 
tures have justified the many studies performed with the unique 
capabilities of Spitzer/IRS. We have used this experience to fo- 
cus one of the last observations of "cold" Spitzer to obtain high- 
quality mid-IR spectra of z ~ 2 starbursts. Our sample was 
5/?/?zer-selected in a field with very rich multiwavelength data. 
We aimed to select starburst-dominated ultra-luminous IR galax- 
ies at z ~ 2 based on the presence of a rest-frame 1.6 pm stel- 
lar bump- caused by the photospheric opacity minimum in gi- 
ant or supergiant stars (e.g. |John 1988] |Simpson & E isenhardt 



1999; Weedm an et al. 2006) - in the 5.8 pm IRAC band, together 
with very strong PAH emission implied by a high 24 pm flux 
density. Using these "5.8/um-peakers," we sought to address in 
the best conditions the various goals allowed by such studies: 
determination of redshifts from PAH features, accurate enough 
for allowing CO line observations with current millimeter facil- 
ities; detailed examination of the properties of PAH features in 
such high-z sources, especially their luminosity ratios, and com- 
parison to local systems; estimation of the ratio between PAH 
bands and the weak underlying continuum that may constrain 
the modest AGN contribution; and accurate assessment of the 
correlation of PAH luminosities with total and far-infrared lumi- 
nosities (and with the radio luminosity). Silicate absorption at 
9.7 pm may provide additional information about the distribu- 
tion of interstellar dust. 

In this paper, we present the observations and the results 
from IRS spectroscopy of a complete sample of 16 sources se- 
lected by their IRAC and MIPS fluxes. The sample selection and 
the observations are described in Section[2] The method adopted 



3] The re- 
and|5] In 



to reduce these observations is presented in Section 
suits and their analysis are discussed in Sections [4 
Section [6] we discuss the discrimination between starburst- and 
AGN-dominated sources, and we make comparisons with the 
other samples of SMGs. Throughout the paper, we adopt a stan- 
dard flat cosmology with Hq = 71 kms -1 Mpc~', Qm = 0.27, 
and Q A = 0.73 ( |Spergel et al. 2003[ ). 



2. Sample selection and observations 

2.1. Sample selection 

Our sample is based on the 5p/fzer/SWIRE "5.8 yum-peaker" 
sample of |Fiolet et al. (20 09) (hereafter F09). Using the 
2007 SWIRE internal catalogue, we have selected the sixteen 
24 /im-brightest sources of this sample with S > 0.5 mJy. 
These sources obey the "5.8/um-peaker" criteria with S3.(,/j m < 
^4.5^m < > Some of our sources have no detec- 

tion in the 8.0 pm band, but we assume that these sources are also 
"5.8 //m-peakers" if their fluxes at 5.8 pm are greater than the de- 
tection limits at 8.0yum (< 40mJy). To remove low-redshift in- 
terlopers, we also require that the sources be optically faint, i.e., 
Ki™, > 23 ( Lonsdale et al. 2006[ ). Our sample is complete for 



these criteria over 0.5 deg 2 . This area includes the "Lockman- 
North" field obser ved with the V LA (|Owen & Morrison 2 008 ; 
|Owen et al. 2009| ) and Herschel ( |01iver et al. 2010| ). For this 
sample from F09, we have 1.2 mm flux densities (Si.2mm) 
obtained with IRAM/MAMBO [Institut de Radioastronomie 
Millimetrique/Max-Planck Millimeter Bolometer] array ( |Kreysa| 
|et al. 1998] ). 62% (10/16) of our sources are detected with 
S/N > 3 at 1.2 mm (Table [T]) and are thus sub millimeter galaxies 
(SMGs) withS^mm > 2mJy «S li2mm ) = 2.72 ± 0.54 mJy). The 
six other sources have an average (Siimm) = 1.01 ± 0.68 mJy 
(see Table [T] and F09). This sample is thus representative of a 
5/?/fzer-selected subclass of powerful z ~ 2 ULIRGs (SMGs and 



OFRGs, Chapman et al. in prep; Magdis et al. 2010), of which 
there are ~ 2000 in all SWIRE fields that peak at 4.5 or 5.8 pm 
and have $24,™ > 500/uJy. 

Nine of our sources detected at 1.2 mm were observed at 
350 pm by Kovacs et al. (2010| . All these sources have been 
detected, and our 16 sources are also detected at 250-350/mi 
by Herschel (Fiolet et al. in prep; Mag dis et al. 2010} . The 
"5.8 //m-peaker" criterium also implies that our sample is a mass 
selected-sample and has z ~ 2. 
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All of our sources benefit from very deep radio data at 1400, 
610, and 324 MHz from the VLA and GMRT ( |Owen & Morrison| 
2008 ; Owen et al. 2009| Owen et al. in prep.), with observed flux 



densities S20_m > 50 //Jy (Table[T]and F09). 

Our sample differs from the other samples of ULIRGs essen- 
tially in its selection criteria and redshift range. Indeed, Farrah 
|et al. (200"8) and Sajina et al. (200 7) have selected brig ht sources 
at 24 pcca.. However, the sample of Sajina et al. (2007) is built to 
select the reddest objects in the [24 /mi - 8 /mi] vs. [24 /mi - R] 
color-color diagram and is composed of PAH-weak and PAH- 
strong sources. In the comparisons that follow, we consider 
only the PAH-strong sources of this sampl e. |Pope et aL~72 008 ), 
|Menendez-Delmestre et al. (2009) , and |Valiante et al. (2007) 
have built samples of SMGs detected at 850 /mi. Nevertheless, 
the sources from Pope et al. (2008 ) selected with S 24mn_>200/£ 



Jy are fainter than our sour ces (£"24^m > 500 /Jy ). |Menendez^ 
|Delmestre et al. (2009) > andJValiante et al. (2007) have also se- 
lected their sample based on 1 .4 GHz flux density. The sample 
of Desai et al. (2009 1 is relatively similar to ours, with sources 
selected by a redshifted 1 .6 pm stellar bump in the IRAC bands 
and S 24 fim > 500 pi y. These sources are also optically faint and 



satisfy the definition of Dust-Obscured Galaxies (DOGs: Dey 
et al. 200"8) . The sample of Huang et al. (2009) is also selected 
with S24/JIH > 500 p Jy and with color-co lor criteria yielding star- 
forming galaxies. Coppin et al. (2010) have built a sample of 
SMGs selected to be AGN-dominated. While our sample has 
redshifts from ~ 1.7 to ~ 2.3, the other samples span larger red- 



shift ranges. |Shi et al. (2009 1 have selected sources at intermedi- 
ate redshift around z ~ 1 . |Farrah et al. (2008| l observed sources 
with z = 1.3 - 2.2 centered at z ~ 1.7 ("4.5 /mi-peakers"). The 
samples of SMG s from|Pope et al. (2008) ; Menendez-Delmestre 



etal. (2009) a nd|Valiante et al (2007| i7and the DOG sample of 
Desa i et al. (2009[Cs pan redshifts from ~ 1 to ~ 3 (0.6 to 3.6 for 
Menendez-Delmestre et al. 2009 ). The sample of star-forming 
galaxies of |Huang et al. (2009} has redshifts between 1.7 and 2.1 



centered at z ~ 1.9. Finally, the sample of Coppin et al. (2010) 
ranges from z ~ 2.45 to ~ 3.4. 

2.2. Observations 

The photometric redshifts of our targets, based on the IRAC 
bands, lie in the range z ~ 1.5 - 2.7 with uncertainties of +0.5 
(see F09). In this redshift range, the low-resolution (R ~ 60 - 
127) Long-Low observing modules of IRS, LL1 (19.5-38.0/mi) 
and LL2 (14.0-21.3 pm), can cover most of the PAH emission 
features from 6.2 to 11.3 pm in the rest frame. Indeed, the actual 
redshifts have all proved to lie in the range 1.7 < z < 2.3 (see 
Sec.[4TT). 

This project (GO50119, PI G. Lagache) was observed for 
a total time of 53 hours (41 hours on-source) during the last 
cold campaign of Spitzer, on January 20-22, 2009. We have di- 
vided our sample into four groups based on 24 pm flux density. 
For each group, the on-target exposure time is summarized in 
Table[2] We chose these exposure times following the guidelines 
from the IRS 1ST report, based on IRS low-resolution observa- 
tions of two faint sources by Teplitz et alQand the Spectroscopy 
"Performance Estimation Tool" (SPEC-PET). We also used our 
own IRS observations of distant galaxies to check the expected 
S/N. The integration times have been fixed to achieve a S/N 
~ 4 - 5 per wavelength element in LL1 and LL2. The obser- 
vations were done in mapping mode to maximize the rejection 
of "rogue" pixels and (immediately after applying the "skydark" 



calibrations) to limit contamination by latent charge from bright 
objects in our "Long" spectra of faint sources. 

3. Data reduction 

The observed spectra were processed with version SI 8.5.0 of 
the Spitzer Science Center pipeline. The pipeline produces the 
basic calibrated data (BCD), which are corrected for different 
effects (e.g., cosmic ray removal, masking of saturated pixels, 
and linearization of the science data signal ramps). 

Ap plying the IDL code IrsLow, developed by |Fadda et al.| 
(2010 1 and available on-lin^] we have extracted the spectra 
from the BCD. IrsLow is a data reduction code built especially 
for low-resolution IRS observations. This code makes a back- 
ground subtraction and extracts the spectra of the target. The 
background is estimated after the user has manually masked the 
regions of interest of the data. In practice, we mask the posi- 
tion of the sources on the spectra to estimate and subtract the 
background. The spectrum is then extracted using a PSF fitting. 
The code also allows the simultaneous extraction of two different 
spectra in the case of two nearby sources. The details of using 
IrsLow and some examples can be found in Fad da et al. (2010} . 
The output of the program is an ASCII file containing the spec- 
trum and its errors. To facilitate our analysis, we make some 
modifications to each of these files. First, we limit the wave- 
lengths to 20.5 pm through 35 pm for the first order (i.e., LL1) 

- the errors on the flux outside this range becoming quite large 

- and to < 20.5 pm for the second order (i.e., LL2). Second, we 
sort the data by wavelength. Finally, we smooth the spectrum in 
IDL. This filter makes a running mean over five points provid- 
ing thus an effective spectral channel width of 0. 10 pm. To avoid 
artificial improvement of S/N, the errors on the spectra are kept 
unchanged by the smooth filter. 

4. Results: Spectra 

The spectra of our 16 sources are shown in Appendix|A](see 
also Figs. [3] and |B). These spectra show high-S/N (> 5<x) de- 
tections of PAH features at 6.2, 7.7, 8.6, and (in many cases) 
11.3/miforz~2(Table[4]). 

Considering the uncertainties of our data, all our sources 
have MIPS 5 2 4^m from the 2007 SWIRE internal catalogue that 
are compatible with the 24 pm flux density {S 24 jU m,spectra) that we 
have directly measured from the spectra after convolution with 
the MIPS 24 pm filter profile (Fig. |A). For 8/16, the difference 
is less than 10%. These flux densities are reported in Table[T] A 
similar result was obtained by Pope et al. (2008 1. 



4.1. PAH redshifts 



http://ssc.spitzer.edu/irs/documents/irs_ultradeep_memo.pdf 



Applying the method described in |Bertincourt et aL~(2 009 ), 
we have determined redshifts from PAH features (zpah)- This 
method is based on the cross-correlation of 21 templates with 
our IRS spectra. These templates are spectra of different kinds 
of sources. Five templates are dominated by PAHs (Smith et al. 
|2007|). 13 a re less PAH-dominated local ULIRGs from |Armu"s 
|et al. (2007 1, including among others Arp220 and Mrk231. The 
remaining three are a radio galaxy, a quasar, and a Wolf-Rayet 
galaxy. The best estimate of the redshift, zpah, is the average 
of redshifts found for all templates after rejection in cases of 
very high^ 2 Or 2 /SNR> 8) or too low SNR (< 1). These red- 
shifts are reported in Table [3] The estimated rms error of zpah 

2 http://web.ipac.caltech.edu/staff/fadda/IrsLow 
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Fig. 1. Ratio of S 24/jm,spectra deduced from our spectra using the 
MIPS 24 /mi filter, to S 2 4^m from the 2007 SWIRE internal cat- 
alogue, as a function of S24//m- The black square is the mean for 
our full sample. The error bars correspond to 1 cr. 



Fig. 2. PAH spectroscopic redshift z PA h as a function of pho- 
tometric redshift from F09. The dashed lines correspond to 
Alog(l + z) = +0.06 or +15%. 



is on average <Az PAH > = 0.010 (median = 0.007; min = 0.004; 
max = 0.064). We also report in Table[3]the redshift, Zbest> of me 
best-matching template. The small differences between z PA h and 
the redshift for the best-matching template are very consistent 
with the estimated rms error (Az PAH ). We use z PA h in the follow- 
ing analysis because this determination is less biased by tem- 
plate selection than the redshift from the best-matching template 
(<Zbest> = 2.011; median = 1.994; scatter = 0.023). We note that 
Zpah and Zbest are very similar for all of our sample, and that 
our sources are best matched by a specific "PAH" template from 
Smith et al. (2007|>. The average of z PAH is 2.017 ± 0.38 (me- 



dian = 1.997). The IRS redshifts of our sample span a narrower 
range, from 1.7 to 2.3, than the photometric redshifts maybe 
due to the strong S^^m- This result is compatible with the fact 
that "5.8/im-peakers" have a high probability of being domi- 
nated by a strong starburst at z ~ 2 because the stellar bump 
at 1.6 /mi rest-frame is redshifted to around 5.8 //m at redshift 
~ 2 (see, e.g., Lonsdal e et al. 2009"] l. Compared to the photo- 
metric redshifts reported in F09 found using the Hyper-z code 
(Bolzonella et al. 2000), the difference is fairly significant (see 
Fig. |2J. Nevertheless, except for four sources, the photometric 
redshifts differ from z PA h by less than 15%, the boundary for 
catastrophic outliers defined by Rowan-Robinson et al. (2008). 



4.2. Spectral decomposition 

We perform spectral decompositions using the PAHFIT code 
developed by |Smith et al. (2007| l. PAHFIT adjusts some PAH 
features and spectral fines as a function of rest wavelength 
(Fig. B). The PAH features are represented by Dmde profiles 
and the atomic and molecular spectral lines by Gaussians. The 
decomposition also tries to fit a stellar continuum component 
(Ti, = 5000 K) and dust components for different temperatures 
(35, 50, 135, 200, 400, 600, 800, 1000, and 1200 K) represented 



by modified blackbodies. Only the intermediate dust tempera- 
tures make significant contributions to the mid-IR continuum. 
The dust extinction, mostly due to silicate absorption, is also 
taken into account. In our decompositions, we have adopted a 
screen geometry, corresponding to a uniform mask of dust ab- 



sorbing the emission, to model the silicate absorption (Smith 
et al. 2007 1. Figure [3] presents an example of the result of the 



spectral decomposition. The PAH complexes at 6.2, 7.7, 8.6, and 
1 1.3 /mi are clearly visible. As proposed by Smit h et al. (2007} , 



the 7.7, 8.6, and 1 1.3 /mi features are better fit by combinations 
of several features at 7.42, 7.60, and 7.85 //m for the 7.7 //m fea- 
ture, at 8.33 and 8.61 /mi for the 8.6 /mi feature, and at 1 1.23 and 
1 1.33 /mi for the 1 1.3 /mi feature. The results of our decompo- 
sitions are somewhat uncertain, due to the difficulty of estimat- 
ing the silicate absorption. However, this absorption is known 
to be weak to moderate in such starburst sources (e.g. 



et al. 2007 


Farrah et al. 2008 


IWeedman & Houck 2008| |Li & 


Draine 200 




Yan et al. 2005 


i. Our targets appear to have es- 



pecially weak absorption making them a good sample for ana- 
lyzing with this decomposition approach. The spectral decom- 
positions for the entire sample are reported in Appendix B for 
individual sources, and in Fig. |4]for the stacked spectrum of all 
16 sources (see also Fig. [12} , It is seen from individual spec- 
tra and the stack that the level of the continuum between 5 and 
1 1 //m remains modest compared to the PAH features (except for 
L10, Sec.[6TT|). 

In the stacked spectrum (Fig. |4}, a weak line is visible at 
~9.7 /mi and identified as the 0-0 S(3) pure rotational transition 
of molecular hydrogen at 9.665 /mi (J = 5 - 3) in the PAHFIT 
spectral decomposition. This line is also visible in the partial 
stacks of Fig. 12 and may be present in the spectral decompo- 
sitions of 13/16 sources (Fig. [B]). The detection of this line is 
confirmed by the detailed data analysis that we perform in a 
separate paper (Fiolet et al. in prep.). The very large luminosity 
of this line, stronger than the CO luminosity, probably implies 
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Ll /: (Jbserved waveleqth (urn) 
20 25 30 




7 8 9 

Wavelength (um) 

Fig. 3. Spectral decomposition for source L17. The solid green 
line is the fitted model. Blue lines above the continuum are PAH 
features; violet lines are spectral lines. Thin grey lines represent 
thermal dust continuum components. The thick red line shows 
the total continuum (stars+dust). The dotted black line shows 
the extinction (e~ T \ = 1 for no extinction). In this example, 
the top axis is labelled in observed wavelength. The redshift is 
Zpah = 1.959. Spectral decompositions for the other sources are 
presented in Appendix B. 



very strong shocks to excite its upper level ~2500K above 
the ground state of H2, although H2 excitation could also be 
influenced by X-rays or UV from AGN. We also see hints of 
other lines, although not at significant levels: the S(5) and S(4) 
transitions of molecular hydrogen at 6.86 /mi and 8.08 /mi, 
respectively, as well as the ionized species [Arlll] at 8.98 /mi 
and [SIV] at 10.52 /mi. 



Stack of all sources 



e 0.2 




~s 0.1 ~ 



Wavelength (um) 

Fig. 4. Spectral decomposition for the stacked spectrum of all 16 
sources. The solid green line is the fitted model. The blue lines 
above the continuum are the PAH features. The narrow violet 
lines are the spectral lines. The grey lines represent the thermal 
dust continuum components. The thick red line shows the total 
continuum (stars+dust). The dotted black line shows the extinc- 
tion (e~ Tj , = 1 if no extinction). 



For each spectrum, we evaluate the continuum flux den- 
sity and luminosity density at 5.8 //m (vL v (5.8 //m)). These val- 
ues are reported in Table [5] There is only a small dispersion 
of vL v (5.8/mi) for our sources (average = 1.87 x 10 11 L Q , me- 
dian = 1.77xlO n Lq, min = 0.71X10 11 L , max = 2.99xlO n L ). 
In Fig. [5] we plot PAH luminosity at 7.7 /mi (Ljj^ m ) as a func- 
tion of vL v (5.8 11m). The PAH luminosity at 7.7 /mi is calculated 
as the integral of the total flux between 7.3 /mi and 7.9 /im af- 
ter subtraction of the continuum contribution (see Eq. [TJ. We 
have also plotted the sample of |Sajina et al. (2007| l. To make a 
true comparison, we have applied a similar reduction and spec- 



tral decomposition to the PAH-rich sources from |Sajina et al. 
(2007|. The values of vL v (5.8/mi) found in this case show no 
significant difference from the published continuum luminosi- 
ties. We see, in Fig. [5] that our sample shows somewhat weaker 
vL v (5.8 /mi) than the PAH-rich sample of |Sajina et al. (2007} or 
the AGN-dominated sample of Coppi rTet al. (2010| l, while the 
values of Lnriym stay in the same range. This is likely due to a 
stronger AGN contribution in the sample of Saji na et al. (20Q7) > 



Note also that the ratio L7.7 ^m/v L v (5.8 /mi) is very consistent 
with that found for the local starbursts of Br andl et al. ( 2006). 
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Fig. 5. Lnnjisa. as a function of the continuum luminosity 
vL v (5.8/mi). The red diamonds are the sources of our sam- 
ple. The open (full) black triangle is the stack of our sources 
with EW7.7^m > (<)6/mi. The full black circle is the com- 
posite spectrum of the SMG sample from Menendez-Delmestre 
et al. (2009) . The black crosses are the PAH-rich sources from 



Sajina et al. (2007 1. The black asterisks are the AGN-dominated 
sources from |Coppin et al. (2010| >. The black line is the best fit 
of all these samples. The dashed black line is the best fit to the 
local starbursts from Brandl et al. (2006) (black squares). 



4.3. PAH features 

Thanks to the spectral decomposition, we have determined 
the fluxes (F) and equivalent widths (EW) for the PAH com- 
plexes at 6.2, 7.7, 8.6, and 11.3 /mi (Table HJ. The flux is cal- 
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Fig. 6. L7.7 Mm as a function of L^^m- The red diamonds are the 
sources of our sample. The open (full) black triangle is the stack 
of our subsample with EW7 7^ > (<) 6 fim. The full black cir- 
cle is the stacked spectrum of the SMG sample from Men endez^] 
Delmestre et al. (2009| l. The black circles are the intermediate- 
redshift star-forming galaxies fro m|Shi et al. (20 09 . The black 
stars are the sample of SMGs from |Pope et al. (2008 1. The black 
crosses are the PAH-rich sources from Saji na et al. (2007) . The 



black asterisks are the AGN-dominated sources from Coppin 



et al. (2010 '. The full black squares are the sources from Farrah 



et al. (2008 1. The open black squares are the local star-forming 
galaxies from Brandl et al. (2006). The solid line is our best fit to 
all these samples. The dashed line is the best fit found by |Pope| 
|et al. (20081 for me i r SMGs and the local starbursts. 

culated by integrating the spectra between 6.2 and 6.3 ^m, 7.3 
and 1.9 fim, 8.3 and 8.7 fim, and 11.2 and 11.4/mi, respectively, 
after subtracting the continuum flux. We calculate the associ- 
ated luminosities L 6 , 2lim , L 7Jlim , L 8 .6//m, and Lii. 3 ^ m , reported in 



Table [5] using: 



AnDlFi 



(D 



where D L is the luminosity distance and F, = f ~(S ^ - S c )dA 
with S ,\ and S c the total and the continuum flux densities at the 
considered wavelengths. 

The equivalent width is calculated following: 



EW 



J,l, 



-dA 



(2) 



All our sources show L 11)lm ^ 5 x 10 10 L and EW 7 7/ml > 2 fim, 
i.e., well above the EWjj m > 0.8 jum limit for "strong" PAH 
sources defined by Sajina et al. (2007 ). All our sources except 
L10 also show EW 62 ^ m » 0.2 //m and L(, 2lim » 2x 10 10 L o . 
These values are usually used to separate starburst- from AGN- 
dominated sources ( Arm uTet al. 2007) . 

Our sample has PAH luminosities on average comparable to 
those of the PAH-rich sources of Sajina et al. (2007 1 and slightly 
larger than those of the SMGs observed by |Pope et al. (2 008 ) 



Fig. 7. L7.7 ^ m as a function of Ls.e,™- The symbols and the line 
are the same as in Fig. [6] 
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Fig. 8. Lj.jpm as a function of Lnj^m. The symbols and the lines 
are the same as in Fig. [6] 



and [Menendez-Delmest re et al. (2009) 1, and of the "4.5 /jm- 
peakers" of |Farrah et al. (2008| l. SucTTa behavior might be nat- 
urally explained by stro nger £24^01 than SM G samples samples 
and higher redshift than |Farrah et al. (2 008 )'s sample. The S/N 
of our average luminosities are higher than the other samples 
(Table [6ji. Indeed, the mean values of L^ ^^ Lj.i/jm, 
and L n . 3jUm are 2.90 ± 0.31, 10.38 ± 1.09, 3.62 + 0.27, and 
2.29 ± 0.26 x 10 10 L o , respectively, for our sample. The stated 
uncertainties are the standard deviations of the means. These val- 
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ues of luminosities are simple averages. The weighted averages 
and the measurements from the stacked spectra of our sources 
give luminosities consistent with the simple averages within the 
uncertainties. For example, the weighted average of L^^jim is 
(2.70 ± 0.09) x 10 10 L , and the value derived from the stacked 
spectrum is (2.81 ± 0.94) x 10 10 L o . 

As seen in Figures [6] [7] and [8] and in Table [7] our PAH lu- 
minosities L6.2/jm, ^7.7//m> ^8.6/an> and yxa have ratios similar 
to those of SMGs, as well as local starbursts and star-forming 
galaxies at z ~ 1. The AGN-dominated sample from Coppin 



et al. (2010 1 seems to have slightly different ratios, with a higher 



The samples plotted in Figs.[6}|7] and[8]follow the relations 

log(L 7 . 7/ , m ) = ( 1. 04 + 0.03) log (L 6 . 2 „ m ) + (0.10 ±0.32) (3) 
log (L 7 . 7ftm ) = (0.87 ± 0.04) log (L 8 . 6/Jm ) + (1.91 ± 0.34) (4) 
log(L 7 . 7/im ) = (1.06 ±0.03) log (Aupm)- (0.03 ±0.30) (5) 

Note that these relations are just plain fits of the displayed sam- 
ples of Figs. [6]-[8] As these samples may be biased in various 
ways, and the properties of the galaxies may vary with redshift, 
one should not expect a strong statistical consistency for these 
relations. The relation s for £6.2 and are consistent with 

the relations found by |Pope et al. (200"8] l for local starbursts and 
SMGs. 



These relations 



-'1.1 pm 



/Ln3fim are 



show that both 
4 for our sources 



Luftm/L^/im and 
at z ~ 2 (see also 

Figs. 13 and 14 and Ta ble |7}, i.e., about twi ce as large as for 
the local starbursts of |Brandl et al. (2006) and the sources 
of |Farrah et al. (2008] > at z ~ 1.7. On the other hand, the 



ratio Isjjn m /Lg£n m is about twice as small in these luminous 
galaxies at z ~ 2 as in fainter local starbursts (Tabl e [7). This 
behavior is similar to that of SMGs (Table [7] and Sec. |6.2) >. The 
difference with respect to local starbursts could be explained by 
a difference in the size or ionization of PAH grains (Sec. |6.1 



|6.2| and Fig. [14}, or by a different extinction due to silicates 
around 10//m (affecting the 8.6 and 11.3 fim bands) or ice at 
~6 fim (absorbing the 6.2 fim band) (Menendez-Delmestre et al. 
2009; |Spoon et al. 2002) . However, the effects of extinction 
seem somewhat contradictory for the different band ratios. 
The finding of a lower 7.7/8.6 ratio compared to normal, local 
starbursts, is consistent with a low silicate extinction in our 
sample, as does the generally low level of extinction in the 
spectral fits. However, the 'larger' 7.7/11.3 ratio that we report 
is going in the other direction; similarly, our 'larger' 7.7/6.2 
ratio would imply a larger ice absorption than in local starbursts. 

L10 is the only source that seems not to follow the rela- 
tion between L 7 7/jm and L6.2//m, showing too weak a feature at 
6.2 fim. However, Table [5] and Fig. [6] show that the deviation of 
L10 from the correlation remains within the statistical errors of 
this noisy spectrum. The offset is probably due to higher extinc- 
tion or a more energetically important AGN (Rigopoulou et al 
[1999) ; see also Sec.|6TT] 

4.4. Silicate absorption 

The determination of the optical depth of the silicate ab- 
sorption at 9.7 fim (T9/7), reported in Table |4] is a product of 
the general fit of the spectra with PAHFIT. It is based on inter- 
stellar extinction similar to that of the Milky Way ( Smith et al. 
|2007| . We note that the silicate absorption is weak on average, 
(rg 7 ) = 0.28 (0.50 if we consider only the sources with Tg 7 + 0) 
vs. 1.09 for the PAH-rich sample of |Sajina et al. (2007[ ). Indeed, 
for 44% (7/16) of our sources, PAHFIT provides t 9 . 7 = 0. It is 



clear that all values of extinction (rg 7 ) are very uncertain be- 
cause of the weakness of the continuum and its uncertainty, and 
because of the lack of measurements at long wavelengths. We 
thus state only the qualitative conclusion that the average extinc- 
tion is likely weak in most of our sources. Considering the values 
of EW6.2 and T9.-7, most of our sources belong to class 1C of the 
diagram of Spoon et al. (2007 , their Figs. 1 & 3). This class is de- 
fined by spectra dominated by PAH emission with weak silicate 
absorption. 



5. PAH emission and star formation 

As noted above, it is well known that infrared PAH emission 
is strongly correlated with star formation through excitation and 
IR fluorescence of PAHs induced by the UV radiation of young 
stars. The exceptional quality of our PAH spectra at z ~ 2 allows 
us to further assess the correlation between PAH luminosity and 
total IR luminosity, and to compare with other starburst samples 
at z ~ 2 and at lower redshift. The deep radio data allow a precise 
comparison with the radio luminosity which is also known to 
trace the starburst IR luminosity well. 

5.1. The correlation between IR and PAH luminosities 




Fig. 9. Ljj^ m as a function of Lir. The symbols are the same as 
in Fig. [6] The solid line is our best fit to all samples. The dashed 
line is the best fit found by Pope et al. (2008 1 from fitting both 



their SMGs and the local starbursts. The dotted line is the best 



fit found by|Shi et al. (2009|>. The dot-dash ed line is the best fit 
found by Menendez-Delmestre et al. (2009). 



The correlation between the strength of the PAH emission 
and the far-IR luminosity in starbursts is well established al- 
though not understood in detail. Remarkably, it has been verified 
among local starburst galaxi es over two orders o f magnitude in 
Lir, from 10 10 to 10 12 Lo, by |Brandl et al. (2006) . The approxi- 
mate proportionality between L 7 . 7/im and Lir has been extended 
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through the ULIRG range up to 10 13 L Q and z > 2 by several 
studies (|Shi et al. 2009, Pope et al. 2008| |Menendez-Delmes"tre 



et al. 2009| l. Our sample gives us an opportunity to check how 
tight the relation is among the strongest starbursts at z ~ 2. 

For the nine sources of our sample observed at 350/zm by 
Kovacs et al. (2010| i - all SMGs (F09) - we have an accurate 
estimate of the infrared luminosity L JR computed by |Kovacs| 
et al. (20101, who fit a distribution of grey -body models with 



different dust temperatures (multi-temperature model). All these 
sources are powerful ULIRGs with L lR close to 10 13 Lq_ (aver- 
age = 8.87 x 10 12 Lq, <AL ir ) = 2.81 x 10 12 L ; see Table|5}. The 
average value of the ratio of IR and PAH luminosities is 



(Ljj^/Ljs) = 11.7 X 10" 



(6) 



It is seen in Table 17] and Fig. [9] that the L 77jUm /LiR ratio is 
within a factor ~ 2 of those for the other high-z samples. 
However, it is a factor ~ 3.5 larger than for less luminous lo- 
cal starbursts as yielded by Eq. [6] Schweitzer et al. (2006 ) have 
also demonstrated the existence of a correlation between the 
strength of the PAH emission and the far-IR luminosity for local 
starburst-dominated ULIRGs and QSO hosts. Nevertheless, the 
(L7.7 tan/Ljs) ratio that they found is probably lower because of 
the more compact nature of those objects. 

In Fig. [9] we plot L 77//m as a function of Lir for different 
samples. We can see that our sample has Lir comparable to 
those of the other SMG samples (Pope et al. 2008; Menendez- 
Delmestre et al. 2009| l, and the PAH-rich sources from Sajina 
et al. (20071 1. The star-ifor ming galaxies of |Shi et al. (2009| > have 
Lir values slightly lower than our sample. We also note that the 
best fit of all these samples is consistent with the relations found 
by |Pope et al. (2008| l and |Menendez-Delmestre et al. (2009 1 
for SMGs and local starbursts (|Brandl et al. 2006) 1 . The AGN- 
dominated sample of |Coppin et al. (20 1 Q] > shows Lir and L 77/im 
greater than our sources, but their ratio is not very different from 
our sample. The sources of Shi et al. (2009 1 seem to follow a 
slightly different relation but with the same slope as the fit found 
by Menendez-Delmestre et al. (2009). The best fit of all these 
samples, from local starbursts to SMGs, is 



log (L 7 . 7 „ m ) = (1.06 ± 0.04) log (Lir) - (2.57 ± 0.52) 



(7) 



This relation is consistent with L 7 7/lm being a good tracer of 
high-z starbursts: a strong L 7 7|Um should be the sign of a high 
SFR. If we use our relation between L 77/im and Lir and assume 
the relation from Kennicutt (1998) between SFR and L JR based 
on a 0. 1-100 M Salpeter-like IMF, we obtain the following re- 
lation 



log (SFR) = 0.94 log (L 7 . 7/im ) -7.27 



(8) 



where SFR is in units of M Q yr 1 . 

This relation is slightly different from that found by 
Menendez-Delmest re et al. (2009| l. The latter is based on the 
luminosity at 7.7 yum calculated from the total flux between 
7.1 jjm and 8.3 /im after continuum substraction, instead of be- 
tween 7.3 /im and 7.9 yum as in our case. However, the values of 
Ljjfim computed with the two definitions are very comparable. 
Applying Eq.[8] we find on average (SFR) = 1260+470 M Q yr -1 . 
This value is consistent with the average SFR of 1060M Q yr _1 
computed from the radio luminosity in F09, and close to the 
median SFR of ~ 1200 M B yr' found for SMGs by |Menendez^ 
IDelmestre et al. (2009| l. 



5.2. The correlation between 1.4 GHz and PAH luminosities 

The correlation between radio luminosity (Li .4ghz) and far- 
infrared luminosity in star-forming regions and in local star- 
bursts is well known (e.g., |Helou et al. 1985[ [Condon 1992; 
ICrawford et al. 1996[ |Sanders & Mirabel 1996|). This cor rela- 
tion has been confirmed at z ~ 2 by |Ch apman et al. (2005 | (see 
alsoF09; ICondon 1992t|Smail et al 2002[|Y un & C arilli 2002t 
Ivison et al. 2002[|Kovacs et al. 2006||Sajina e t al. 2008 1 | Younger] 



et al. 2009[|Ivison et al. 2010[|MagnelIi et al. 2010 1. 
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Fig. 10. Ljjfim as a function of rest-frame vL v (1.4GHz). The 
symbols are the same as in Fig. [6] The black line is our best fit 
to all these samples. 



From the correlation between Lir and L 7 7A(m (Fig. [9] Eq.|JJi, 
we might expect a very good correlation between PAH and radio 
luminosities, as previously observed by Po pe et al. (2008[ ). Since 
our sources benefit from exceptionally good radio data at 20 cm, 
50 cm, and 90 cm (Sec. 



2.1 



it is worthwhile to directly check 
the tightness of this relation. 

It is straightforward to compute the rest-frame luminosity at 
1.4 GHz (Li.4ghz = vL y (1.4GHz)) using, e.g., the formula given 
by |Kovacs et al. (2006) : 



vL v (1.4GHz) =47rL>£s 2 0cm(1.4x 10 9 )(1 +z) 



(-ff-i) 



(9) 



where D L is the luminosity distance and a is the radio spec- 
tral index defined by the best power law fit, S v oc y", between 
520cm, 550cm, and Sever n- We have computed vL v (1.4GHz) fo r 
our sample and those of Sajina et al. (2007 ), Pope et al. (2008), 
Shi et al. (2009] ), and |Brandl et al. (2006j l For our sample, we 
have used a as computed in F09 and reported in Table [T] For the 
PAH-rich sample of Sajina et al. (2007 ), we have used the val- 
ues published in Sajina et al. (2008[ l. For the other samples, we 
use the average value of a, (a) = -0.64 found in F09, which is 
consistent with the typical radio spectral index for star-forming 
galaxies (e.g., |Condon 1992 1. The values of vL v (1.4GHz) are 
reported in Table [5] 
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As expected, we observe a correlation between L\ .4 GHz and 
Lin nm- Considering our sample, but also the local star-forming 
galaxy sample (Br andl et al. 2006) , the sample of SMGs (|Pope| 
|et al. 2008) , the PAH-rich sources from Sajina et al. (2007) , and 
the intermediate-redshift starbursts ( jShi et al. 2009| l, we obtain 
the following relation (Fig 10 1: 



log(L 7 . 7/im ) = (0.91+ 0.06) log (L v ( 1 .4 GHz)) + (4.60 ± 0.36)( 1 0) 

The effect of the positive ^-correction and the value of a are par- 
ticularly important for the determination of Li.4ghz and ma y sig- 
nificantly affect the reliability of the relation stated in Eq. 10 The 



correlation between Li.4ghz and L 77/im seems to be more scat- 
tered than the correlation between L IR and L 7 . 7/Jm . It is not sur- 
prising considering the uncertainties on Li.4ghz and the scatter of 



the correla tion between L\ R and L\ 4gh z as shown by e.g. |Kovacs 
|et al. (2010| i. 

5.3. PAH emission and stellar mass 
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vL T (1.6^tm) [L e J 



Fig. 11. L-irifmi as a function of the stellar bump luminosity 
Li.6/mi- The symbols are the same as in Fig. [6] Only the sam- 
ples at z ~ 2 are plotted. The line is the bestnt found for the 
"normal" star-forming galaxies of |Lu et al. (2Q03) >. 



The stellar bump luminosity at 1.6 pm (vL v (1.6/im), re 
ported in Table [5] gives an estimate of the stellar mass (M*) 



of our sources (see, e.g., F09; Seymour et al. 2007 Lonsdale 
et al. 2009) . Applying the same mass-to-light ratio as F09 



(M + /vL v (1.6jum) = 0.7 M Q /L ), M* derived from vL v (1.6pm) 
ranges from 1.1 to 3.0 x 10 11 M G (average = 1.6xlO n M Q , me- 
dian = 1. 4xl_0' 1 M Q ). 
In Fig. 
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we plot L 77jUm as a function of v L v (1.6 pm) for 
our sources and the PAH-rich sample from Saji na et al. (2007). 
Our sample shows no clear difference from the sources of[Sajina 
et al. (2007[) in terms of v L v (1.6pm). Both samples follow 



the relation found for "normal" star-forming galaxies (Sajina 



vL v (\.6pm) do not trace AGN (Saji na et al. 2007) . Since the 
PAH luminosity is tightly correlated with the star formation rate, 
the relation between vL y (l.6pm) and L 77/im seems to confirm 
the existence of a correlation between star formation rate and 
stellar mass at z ~ 2 as previously observed by Pannella et al. 
|(2009) and |Mobasher et al. (200"9) . However, note that the mass- 
to-light ratio remains uncertain because of uncertainty in the age 
of the dominant stellar population. 



6. Discussion 

6.1. Distinction between starburst and AGN 
Stack of sources with EW>6um 



E 0.2 




6 8 10 

Wavelength (urn) 

Stack of sources with EW<6um 



E 0.2 r 




et al. 2007 Lu et al. 2003 i. This result implies that L 77/Jm and 



Wavelength (urn) 

Fig. 12. PAHFIT spectral decomposition of stacked spectra. 
Upper panel: 7 sources with EW 7 7A , m > 6/im. Lower panel: 9 
sources with EW 7 7//m < 6 pm. The solid green line is the fitted 
model. The blue lines above the continuum are the PAH features. 
The narrow violet lines are the spectral lines. The grey lines rep- 
resent the thermal dust continuum components. The thick red 
line shows the total continuum (stars+dust). The dotted black 
line shows the extinction (<?~ Tl = 1 if no extinction). 



As previously discussed, the relative strength of the PAH 
features and the mid-IR continuum is a very efficient crite- 
rion for distinguishing starburst- from AGN-dominated ULIRGs 
(e.g., |Yan et al. 2007||Sajina et al. 2007] |Weedman et al. 2006[ 
Houck et al. 2005) . |Genzel et al. (1998) and |Rigopoulou et al. 
(1999) use the parameter 7.7 pm-(L/C), defined as the ratio be- 
tween the peak flux density at 7.7 pm and the continuum flux 
density at the same wavelength, to distinguish between AGN 
and starburst. AGN-dominated composite sources show values 
of 7.7/mi-(L/C) weaker than starburst-dominated systems. An 
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equivalent straightforward indicator is the equivalent width of 
the PAH features, especially EWjj^ m , since it is proportional 
to the ratio of L-u^ to the 7.7 //m continuum. The equivalent 
widths of our sources span a broad range, more than a factor of 5. 
Despite the large uncertainty in the continuum, e.g., as estimated 
by PAHFIT, it is clear that EW7 .7^,^ should be a relevant param- 
eter for discussion of AGN contributions. For our entire sample, 
we find (EW^v^m) = 5.64 /im (median = 4.88) vs. 3.04 yum for 



the stacked spectrum of Menendez-Delmestre et al. (2009 ). 

Based on the values of EWyj^m (see Table |4J, we have 
constructed two subsamples: the first has EWjj^m > 6 fim (7 
sources, (EW^^) = 8.33 /im) and the second has EWjj^ < 
6 fim (9 sources, (EWjj^ m ) = 3.55 fim). We have produced 
stacked spectra of these two subsamples, built from the weighted 
means of the individual spectra scaled to (Fj j Mm ). They are plot- 



ted in Fig. 12 It is obvious, as expected, that the second subsam- 



ple has a significantly stronger continuum. 



Applying the same method as Menendez-Delmestre et al. 
(2009), we estimate an X-ray luminosity (L x ) from the 10 fim 



continuum flux (Table [4]l via fhejsKrabbe et al. (200 l ) relation. 
Then assuming the Lx/Z-ir found for AGN by Alex ander et al.| 
(2005), we estimate the contribution of AGN to Lir. We find 
that the AGN contribution is ~ 20% for our whole sample and 
~ 10% and ~ 25% for the subsampl es with EWjj^m > 6 fim and 
with EWjj Mm < 6 fim, respectively. [Menendez-Delmestre et al.l 
|(2009| ) have found ~ 32% to be the AGN contribution in SMGs. 
The difference can be interpreted as a difference in the contin- 
uum flux density, with continuum larger by a factor ~ 2 for the 
sample of Menendez-Delmestre et al. (2009). 

In Fig. [9] (see also Table |6J, it is seen that the two subsam- 
ples are not different in term of L\r and follow the correlation 
between Lj j^ m and L\r found in Sec. 5.1 This result tends to ar- 
gue that the majority of the IR luminosity is associated with the 
starburst and not with the AGN. This is not the case in the sam- 
ple of Menendez-Delme stre et al. (2009| ), which shows a smaller 
Lu^m/LiR ratio, implying that the AGN contributes more to Lir 
(Table|7]). 

Another possible sign of a strong AGN contribution is a 
large Lin ^ mlL^iym ratio, as di scussed by, e.g., |Rigopoulou et al. 
( 1 999j) and [Pope et al. (2008[ ). They have observed such a rise 



in Ljj Mm versus Z-6.2/jm for some SMGs, which they attribute to 
greater silicate absorption in the presence of a more luminous 
AGN than is seen in most SMGs or starburst galaxies. Stronger 
silicate absorption significantly affects Ls.e/jm and Lhj^, while 
it may lead to an overestimate of Ljj^. However, as seen in 
Figs. [6] and [8] there is no significant difference in the value of 
£7.7 /jml 'L(,.2ijm between sources with large and small EWs, per- 
haps because the silicate absorption is small. 



We see in Fig. [TO] and Table [6] that sources with EW 



7.7 fim 



6 /im seem to have an average value of vL v (1.4GHz) somewhat 
greater than that of sources with EW7 7^ > 6 fim. Radio ex- 



cess could be a good indication of AGN strength (Seymour et al. 
[2008) [Archibald et al. 20"0Tj |Reuland et al. 2004| ). However, we 
note that practically all of this difference is due to two sources in 
the small-EW subsample, Lll and L17, which are close to the 
radio loudness Unfi t at z ~ 2 as defined by Jiang et al. (2007]) and 
Sajinaetal. (2008|): vL v (1.4GHz) = 3.66x10' L Q . They are thus 



good candidates for having significant AGN emission. Applying 
the same method as previously to determinate the AGN contri- 
bution, we find 26% and 19% for Lll and L17, respectively. 

To summarize: although small, it seems that the AGN con- 
tribution is slightly greater in the sources with EW 77jUm < 
6 fim than the systems with EWjj^ > 6 fim. Nevertheless, all 
these sources show strong PAH emission and individually fol- 



low more or less the relations found for local starbursts ( |Brandl| 
|et al. 2006) . They are probably all starburst-dominated or AGN- 
starburst composite sources, with at most a weak AGN contribu- 
tion to the mid-infrared. 

Our sample is particularly homogeneous in term of PAH 
luminosity (see Table [5]). Only one source, L10, is different. 
Indeed, this source shows a spectrum slightly different from the 
others. Its continuum is stronger and has a larger slope (Fig.[B|. 
Moreover, L10 deviates from the correlation between Z^^m 
and L-irijim with a low L^.i^m- L10 also presents the weakest 
EW^^m, EW 77jUm , and EWg 6/jm of our sample, with values 
of 0.10, 2.32, and 0.77 fim, respectively. Such a difference can 
be explained by a larger AGN contribution to the mid-infrared 
emission, and also by a low signal-to-noise ratio (see Fig [A| . 
Nevertheless, the fact that L10 follows the other correlations - 
especially the correlation between vL v (1.4GHz) and L-jj^ - 
seems to prove that this source is starburst-dominated. The same 
kind of source was observed by |Pope et al. (2 008 ). 

6.2. Comparison with other studies 

The present work is based on a complete sample of 16 
sources selected to be starburst "5.8 /mi-peakers" that are bright 
at 24 fim, with S 24 fim > 0.50 mJy. The average flux density of 
our sample is (S 24 fim) = 0.63 mJy. This is a brighter sample 
of SMGs co mpared to that of |Pope et al. (2008| ) ((S 2 4jUm ) = 
0.45 mJy) or |Menendez-Delmestre et al. (2TO9J «5 2 4 ^m) = 
0.33 mJy). On the other hand, the PAH-rich sample of|Sajina 
|et al. (2007 1 is biased toward the brightest 24 fim sources, with a 
mean (S24fim) = 1.19 mJy. Our sample is also brighter than the 
ULIRGs at z ~ 2 of |Fadda et al. (2010[ ) «5 2 4 /Jm ) = 0.27 mJy). 
All these samples have an average redshift z ~ 2, but ours spans 
a narrower redshift range than the others. Indeed, our sample 
has zpah between 1.75 and 2.28, while the samples of [Pope et al. 
(200 8) and|Menendez-Delmestre et al. (2009| ), the ULIRGs from 



|Fadda et al. (2 010) and the PAH-rich sample of |Sajina et al. 
|(2007| ) span z~^0.93 - 2.55, 0.69-3.62, 1.62-2.44, and 0.82- 
2.47, respectively. 

The differences between PAH luminosities for all the sam- 
ples are small, but our PAH luminosities show an higher S/N 
than the other high-redshift samples (Table [6}. Our sample has 
PAH luminosities in all bands quite comparable to those of 
Sajin a et al. (2007|) and about twice as large as th ose of Pope| 
|et al. (2008| ) and |Menendez-Delmestre et al. (2009| ). 

As regards the PAH-band luminosity ratios, we see in Table[7] 
that our sample is not very different in term of ratios from the 



other SMG s amples (|Pope et al. 2 008 
et al. 2009[ |Sajina et al. 2007|) 



Menendez-Delmestre 
Nevertheless, IMenendez- 



Delmestre et al. (2009) find lower ratios for L^iim/Lm and 
^7.7/W^ir- This is probably due to a greater AGN contribution 
to L IR (Sec.[6~T). 

Our sample also sho ws a higher L l llim IL iR ratio than the 



lower-redshift sample of Brandl et al. (2006). This difference 
could be explained by a more extended PAH distribution in our 
sources. Indeed, if the PAH distribution is extended, a larger 
fraction of PAH can survive in strong UV radiation field from 
AGN or young stars. Then Ljj^/Lir becomes higher than in 
compact PAH distributions where PAH are easier destroyed by 
UV radiation field (e.g., Huang et al. 2009[ ). 

In Fig. [T3] we plot the different PAH ratios of these samples 



following O'Dowd et al. (2009| ). We find that our sample and 

6.2 frni/L j 1 .3 ^ 



the other SMG samples show L^ liim jL 



7.7 /a n /Ml. 3 <nn 



VS. L fi 



ratios consistent with the sample of O'Dowd et al. (2009). The 
sample of Brandl et al. (2006 ) falls in the border of the dispersion 
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for the PAH ratios of Galactic HII regions, dwarf spirals, and 
starburst galaxies found by Galliano et al. (2008 1. The plot of 
Lj.inm/Lu,3tim vs. U. 6lim /L n3fim shows that all these samples 
have ratios different from those of |Q'Dowd et al. (2009[ ) because 
of stronger emission at 8.6jum. Our sample does not follow the 
relation found by Galliano et al. (2008). 
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Fig. 14. Adapted from O'Dowd et al. (2009[ >: L U3/jm /L 7 . llJm as 
a function of L6.2/W 'tan ^m- The symbols are the same as in 



Fig. 13 The black lines represent the expected ratios for fully 
neutral or fully ionized PAH molecules of a given number of 
carbon atoms from the models of Drai ne & Li (2001) . 

A likely explanation of these differences in PAH ratios is 
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/j.m/ 11.3 fim 



modification of the size or ionization of PAH grains. Fig 
adapted from |Q'Dowd et al. (2009[ ), seems to prove that our 
sample has PAH grains with sizes comparable to those in the 
other samples of SMGs, but with a slightly lower fraction of 
ionized grains than the sampl es of |Pope et al. (2008] l and (es- 
pecially) Sajin a et al. (2007] l. |Brandl et al. (2006) find smaller 
and more neutral grains in local starbursts than our sample. The 
difference could arise due to the strongest ionizing UV radia- 
tion from young stars in the most powerful starbursts but less 
ionising than in the sample of Pope et al. (2008 1 and Menendez- 
|Delmestre et al. (20"09] > probably because of lower AGN contam- 
ination. 



Fig. 13. Upper panel: Ljj m /Ln.3pn 
.2^m/^n.3/im- Lower panel: L-j j Ln 3 ^ m 



a function of 
as a function of 



L$.6iim/Lu.3iim- in both panels, the red diamonds are the sources 



in our sample. The small black circles are the sample of O'Dowd 
et al. (2009) . The open (full) black triangle is the stack of our 
sources with EW 77A(m > (<)6/mi. The black cross shows the 
average for the PAH-rich sources from Saji na et al. (2007| ). The 



black star shows the average for the sample of SMGs from Pope 
et al. (2008] l. The open black square shows the average for the 



local star-forming galaxies from Brandl et al. (2006). The full 



black square is the average for the sample of "4.5 jum-peakers" 
from Far rah et al. (2008) . The large black circle is based on 



the composite spectrum of the SMG sample from Menendez 
|Delmestre et al. (2009 ). The black line and the dotted lines 
are the best linear regression and the lcr dispersion found by 
Galliano et al. (2008] l for Galactic HII regions, dwarf spirals, 
and starburst galaxies. 



7. Summary and conclusions 

In this paper, we have presented the results of Spitzer/JRS 
observations of a complete sample of 16 sources selected to 
be 24 jum-bright and "5.8jt/m-peakers." The spectra obtained 
show very strong PAH features at 6.2, 7.7, 8.6, and 11.3/im, 
along with a weak continuum. Thanks to these exceptionally 
strong features, we have estimated accurate PAH redshifts that 
span the range z = 1.75 - 2.28. Our average PAH redshift is 
(zpah) = 2.02 + 0.15, and the average error of the individual zpah 
is (Azpah) = 0.010 (median = 0.007; min = 0.004; max = 0.064). 
The stacked spectrum of our 16 sources displays evidence for the 
pure rotational 0-0 S(3) molecular hydrogen line at A = 9.67 /jm, 
which is confirmed and analyzed in a separate paper. 

Thanks to the very good quality of our IRS spectra, we have 
calculated PAH luminosities L^^jim, J-nniim, i8.6/jm, and Lh j^ 
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to have average values 2.90 + 0.31, 10.38 ± 1.09, 3.62 + 0.27, 
and 2.29 ± 0.26 x 10 10 L Q , respectively. These luminosities have 
a S/N higher than the other SMGs samples. We have studied the 
correlation between Ljj^ m and the other PAH luminosities. All 
our sources, except perhaps L10, follow a correlation similar to 
that found for local starbursts, and they are not very different 
from other samples of SMGs at z ~ 2. We have also confirmed 
the very good correlation, previously observed, between Ljj^ m 
and Lir at high redshift, and we have verified its extension to the 
radio luminosity vL v (1.4GHz). It also extends more loosely to 
stellar mass through vL v (1.6 fim). 

All these relations, the luminosity ratios, and the equiva- 
lent widths allow us to estimate the AGN contribution to the 
mid-IR luminosity. We conclude that our sources are starburst- 
dominated and that the AGN contribution is ~20%. This sample 
is the most pure selection of massive starbursts at z ~ 2 com- 
pared to the other z ~ 2 Spitzer-selected samples. The known 
fact that the equivalent width of PAH features is a good discrim- 
inant between starburst-dominated and AGN-dominated sources 
is reinforced by the study of two subsamples of sources with 
EWvj^m > 6/im and EWjj^ < 6//m. The subsample with 
EW 7 7flm < 6 yum shows a larger AGN contribution (~ 25%) than 
the EW 7 . 7jUm > 6yum subsample (~ 10%), although the AGN 
contribution is small in both subsamples. 
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Appendix A: Spectra and best fit-templates 

In Fig. |Aj we present our spectra and the best-fit templates 
among the 21 templates described in Sec 4.1 using the method 
described in Bertincourt et al. (2009). We also plot the 24 //m 
flux densities from the SWIRE catalogue and the values deter- 
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Fig. A.l. IRS spectra of our sources. The black line is our data. In shaded grey, we plot the l<x deviation. The green short-dashed 
line is the best-fit PAH template from Smith et al. (20071. The red diamond is the 24 /mi flux density from the SWIRE catalogue. 
The red square is the 24 /mi flux density extracted from the spectrum with the 24 /mi filter profile. The pink shape shows the region 
where the transmission of the 24 /mi is maximum. All these spectra are plotted in observed wavelength. The templates are scaled to 
the "(1 + z) x 7.7 /mi" flux. The source ID, the names of the best-fit templates, and zpah are reported in each panel. 
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Fig.A.l. Continued 
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mined from the IRS spectra convolved with the MIPS 24 //m fil- 
ter. 



Appendix B: Spectral decomposition of the spectra 



In Fig.|B] we present the decompositions of our spectra made 
with the PAHFIT code developed by|Smith et al. (2007|>. 



16 



N. Fiolet et al.: Mid-infrared spectroscopy of Spitzer-selected ultra-luminous starbursts at z~ 1 



E 0.2 




0.4 



0.3 



6 


8 10 


12 




Wavelength (urn) 






L9 







E 0.2 r 



0.1 r 



0.0 



-0.1 




7 8 9 

Wavelength (urn) 

L11 



10 11 




7 8 9 10 

Wavelength (um) 




0.0 9 " J ' - Stf"- 



8 10 
Wavelength (um) 




-0.1 h 



7 8 9 
Wavelength (um) 

L10 



10 11 



0.4 



0.3 



E 0.2 r 



-0.1 




0.1 r 

ElnMlfflHitf 'ff Willi! Ill 

ipwrfi rt i 



8 10 
Wavelength (um) 

L14 




7 8 9 

Wavelength (um) 



Fig. B.l. Decompositions of the spectra of our sources. The solid green line is the fitted model. The blue lines above the continuum 
are the PAH features. The violet lines are the spectral lines. The thin grey lines represent the thermal dust continuum components. 
The thick red line shows the total continuum (stars+dust). The dotted black line shows the extinction (e~ T| , = 1 if no extinction). 
The source ID is reported as the title of each panel. 
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Table 1. Spitzer mid-IR, submillimeter, and radio data of the selected sample. 



ID" 


IAU name 


r b 
J 24/im 


O C 

" 24 fan, spectra 


r d 
J 350/jm 


c a 
" 1.2 mm 


S 20cm e 


a' 












0«Jy) 




(mJy) 




0"Jy) 




LI 


SWIRE4. 


J104351. 17+590058.1 


664+ 


18 


674+38 


30.0+10.5 


2.95h 


:0.66 


76.3+9.0 


-0.84 


L3 


SWIRE4 


J104410.00+584056.0 


589+ 


:12 


495: 


t80 




1.04h 


=1.11 


105.5+35.0 


+0.86 


L5 


SWIRE4. 


J104427.52+584309.6 


542+ 


:11 


645: 


t26 




2.75+ 


:0.76 


71.0+20.0 


-1.03 


L7 


SWIRE4 


J104430.60+585518.6 


825+ 


13 


871 + 


:105 




1.39+ 


:0.79 


101.7+13.3 


-0.81 


L9 


SWIRE4. 


J104440.25+585928.3 


653+ 


14 


586+ 


109 


55.1+17.1 


4.00h 


:0.55 


116.5+9.2 


-0.54 


L10 


SWIRE4 


J104549.77+591904.0 


626+ 


12 


605: 


t76 




1.39+ 


:0.81 


67.2+11.7 


-0.40 


Lll 


SWIRE4. 


J104556.90+585318.8 


661+ 


13 


825: 


t55 


38.9+5.8 


3.08+ 


:0.58 


314.8+4.1 


-0.45 


L14 


SWIRE4 


J104638.68+585612.5 


631+ 


13 


661: 


t48 


29.7+4.9 


2.13+ 


:0.71 


159.5+6.0 


-0.73 


L15 


SWIRE4. 


J104656.46+590235.5 


579+ 


:14 


695: 


t41 


25.7+5.9 


2.36+ 


:0.62 


68.6+3.7 


-1.03 


L17 


SWIRE4 


J104704.97+592332.3 


644+14 


629: 


t83 


45.3+12.4 


2.24h 


:0.64 


341.0+33.0 


-0.76 


L19 


SWIRE4. 


J104708.78+590627.2 


556+ 


:14 


618: 


b42 




0.89+ 


:0.97 


102.2+12.0 


-0.44 


L20 


SWIRE4 


J104717.96+590231.7 


692+ 


:15 


558+ 


107 


51.6+7.2 


2.66+ 


:0.78 


51.1+4.7 


-1.25 


L25 


SWIRE4. 


J104738.32+591010.0 


704+ 


13 


667: 


t74 


36.1+7.1 


2.56d 


:0.74 


69.0+9.0 


-1.09 


L27 


SWIRE4 


J104744.60+591413.5 


533+ 


13 


591: 


t50 


27.4+7.9 


2.48+ 


:0.73 


77.2+13.7 


-0.38 


L31 


SWIRE4. 


J104830.70+585659.3 


640+ 


:16 


640+ 


101 




1.62+ 


:0.60 


121.0+39.0 


-0.31 


L33 


SWIRE4 


J104848.23+585059.3 


536+ 


:15 


487: 


t61 




-0.28: 


tl.09 


109.0+33.0 


+0.27 



"IDs and S i 2mm are the same as F09. The bold values are >3 cr detections 
*S , 24 J im comes from the internal SWIRE catalogue. 

C S 24(/m,s P ectra is calcu lated by convolution of the 24 fim MIPS filter and our spectra. 
d S 35o^m comes fro m|Kovacs et al. (2010|>. 
e ^20cm comes from Owen & Morrison (2008 ). 

'The radio spectral index a is calculated using the best power law fit, 5 V oc V, between S20cm> ^50om> an d ^gocm as reported in F09. 



Table 2. Summary of Spitzer IRS observations. 



ID 


LL1 time on tar« 


;et LL2 time on target 






(sec) 


L7 


2x120x6 


8x120x6 


L1,L9, L11,L17,L20, L25 


3x120x6 


8x120x6 


L3, L10, L14, L15, L31 


4x120x6 


9x120x6 


L5, L19, L27,L33 


5x120x6 


11x120x6 



Table 3. Redshifts. 



ID 


Zpah" 


v b 


mid-IR best matching template 
from|Smith et al. (2007) 


LI 


2.200+0.005 


2.196 


pah2 


L3 


1.921+0.005 


1.922 


pah3 


L5 


1.750+0.007 


1.756 


pah3 


L7 


1.944+0.064 


1.886 


pah2 


L9 


2.041+0.007 


2.046 


pah5 


L10 


1.836+0.005 


1.834 


pah3 


Lll 


1.946+0.010 


1.946 


pah5 


L14 


2.032+0.005 


2.034 


pah4 


L15 


1.851+0.010 


1.844 


pah4 


L17 


1.959+0.004 


1.964 


pah3 


L19 


2.240+0.012 


2.222 


pah2 


L20 


2.072+0.010 


2.066 


pahl 


L25 


1.972+0.005 


1.974 


pah3 


L27 


2.022+0.008 


2.014 


pah4 


L31 


2.284+0.009 


2.272 


pah2 


L33 


2.206+0.010 


2.200 


pah3 



"Average PAH spectroscopic redshift. 
6 Redshift for best-matching template. 
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Table 4. Fluxes and equivalent widths for PAH features from PAHFIT ( Smith et al. 2007 1. 



ID 


^6.2 fim 








EW 6 . 2 ,, m 


EW 7 . 7(m EW«. 6fm 


EW n . 3fjm 


^9.7 


S lOiim" 






(10- 22 W.cirr 2 ) 








(yum) 






(mJy) 


LI 


4.93+0.51 


21.69+1.80 


5.48+1.03 




2.86 


9.42 


2.70 




1.08 


0.037 


L3 


3.44±0.51 


8.28+0.69 


3.22+1.30 


2.94+0.85 


2.90 


8.17 


2.67 


1.85 


0.00 


0.036 


L5 


3.82±0.51 


12.41+1.37 


5.16+1.07 


2.48+0.83 


3.67 


12.02 


4.99 


1.63 


0.02 


0.033 


L7 


3.11+0.56 


12.44+1.82 


6.22+2.24 


2.13+1.26 


0.55 


2.82 


1.37 


0.52 


0.00 


0.133 


L9 


2.41±0.38 


8.68+0.67 


4.53+0.61 


1.28+0.91 


0.63 


2.41 


1.33 


0.27 


0.00 


0.096 


L10 


0.29±0.50 


7.99+1.51 


2.83+1.33 


2.10+0.89 


0.10 


2.32 


0.77 


0.48 


0.00 


0.128 


Lll 


3.61±0.48 


12.99+1.74 


6.43+1.70 


3.06+1.68 


0.81 


4.72 


2.28 


1.15 


0.10 


0.080 


L14 


3.37±0.57 


13.17+1.10 


4.39+1.30 


2.50+1.87 


1.19 


7.54 


2.80 


3.17 


0.00 


0.034 


L15 


4.47±0.33 


14.48+1.30 


5.83+0.52 


4.25+0.53 


1.73 


6.88 


3.22 


3.92 


0.00 


0.044 


L17 


2.38+0.51 


14.10+1.82 


5.39+0.75 


3.16+0.96 


0.84 


5.04 


2.21 


1.71 


0.54 


0.051 


L19 


4.93+0.35 


15.67+1.01 


3.87+0.48 




1.71 


7.52 


2.30 




0.00 


0.035 


L20 


3.85+0.63 


9.49+1.49 


3.11+1.77 


5.25+2.46 


1.31 


4.40 


1.30 


2.83 


0.18 


0.066 


L25 


5.39+0.50 


12.45+1.00 


5.98+1.08 


4.82+1.51 


1.72 


2.86 


1.57 


1.43 


1.32 


0.067 


L27 


3.77+0.42 


12.77+1.54 


3.40+0.77 


3.54+1.18 


1.81 


3.43 


1.12 


1.01 


1.09 


0.069 


L31 


4.07+0.52 


16.25+1.37 


4.82+0.60 




1.26 


6.80 


2.44 




0.03 


0.043 


L33 


3.12+0.41 


10.57+2.11 


2.29+0.86 




0.93 


3.93 


1.01 




0.10 


0.054 



"Sioiim is the continuum flux density at lOyum. 
Table 5. Luminosities for different wavelengths. 



ID 


vL v (1.6fim) a 
(10" 


vL„(5.8pm) 


^6.2 fim 


^7.7 jum 

(10 10 




Li 13 pm 


vL v (1.4GHz) a 
(10 6 L o ) 


(1O 12 L ) 


LI 


2.09 


1.55 


4.83+0.50 


21.24+1.76 


5.36+1.01 




8.71+1.00 


11.00+3.80 


L3 


2.01 


0.97 


2.41+0.35 


7.80+5.12 


2.26+0.91 


2.07+0.59 


1.35+0.43 




L5 


1.80 


0.71 


2.12+0.29 


6.90+0.76 


2.87+0.60 


1.38+0.46 


5.75+1.59 




L7 


1.80 


2.99 


2.24+0.40 


9.73+7.18 


4.49+1.62 


1.54+0.91 


8.32+1.15 




L9 


2.00 


2.27 


1.96+0.31 


7.06+0.55 


3.68+0.50 


1.04+0.74 


7.94+0.55 


9.77+3.37 


L10 


1.85 


1.06 


0.18+0.31 


5.01+0.95 


1.77+0.83 


1.32+0.56 


3.16+0.58 




Lll 


1.55 


2.54 


2.61+0.34 


10.77+4.42 


4.65+1.23 


2.22+1.21 


17.78+0.41 


8.13+2.24 


L14 


2.14 


1.96 


2.72+0.46 


11.63+5.30 


3.54+1.05 


2.02+1.51 


13.49+0.62 


8.51+2.54 


L15 


2.97 


1.50 


2.86+0.21 


9.25+0.83 


3.73+0.34 


2.72+0.34 


6.31+0.29 


6.31 + 1.89 


L17 


2.38 


1.57 


1.76+0.38 


10.37+1.34 


3.96+0.55 


2.33+0.71 


26.92+2.48 


9.12+3.36 


L19 


3.63 


2.59 


5.05+0.35 


16.04+1.04 


3.96+0.49 




7.59+0.87 




L20 


4.19 


2.05 


3.25+0.53 


9.82+5.96 


2.63+1.50 


4.43+2.08 


7.94+0.73 


12.30+3.40 


L25 


2.36 


1.54 


4.03+0.37 


9.30+0.75 


4.47+0.81 


3.60+1.13 


7.94+1.10 


7.76+2.50 


L27 


1.74 


1.13 


3.00+0.34 


10.16+1.22 


3.18+0.61 


2.82+0.94 


4.37+0.80 


6.92+2.23 


L31 


1.55 


2.97 


4.37+0.56 


17.44+1.47 


5.18+0.64 




7.94+2.56 




L33 


2.53 


2.43 


3.08+0.40 


10.42+2.08 


2.26+0.85 




3.39+1.01 





a From F09. 

fc From |Kovacs et al. (ZOlot . 
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Table 6. Average luminosities for different samples. 



Sample 



.2 fim 

(mJy) 



^7.7 fim 

(10 10 



-41.3 /jm 



vL„(1.4GHz) 
(1O 6 L ) 



(10 12 L o ) 



This work 


2.017 


0.63 


2.90±0.31 


10.38+1.09 


3.62= 


=0.27 


2.29= 


t0.26 


8.68=t 


1.56 


8.874 


:0.64 


This work, EW 7 7; , m >6/im' 1 


2.040 


0.60 


3.07±0.13 


11.30±0.35 


3.42±0.15 


2.51= 


t0.29 


7.31i 


1.38 


8.614 


= 1.35 


This work, EW 7 .7 U „,<6 fim" 


2.000 


0.65 


1.94±0.14 


7.96±0.72 


2.86= 


=0.21 


2.12+0.29 


9.75=t 


2.59 


9.004 


:0.78 




Sajina et al. (2007}^ 
Menedez-U. et aT(2009) c 


1.524 


1.19 


2.66±0.46 


13.17+3.28 


2.87= 


=0.67 


1.71 = 


t0.37 


18.33+6.50 


6.324 


= 1.36 




2.000 


0.33 


1.09±0.22 


4.10±0.93 


1.63= 


=0.17 


1.39= 


t0.10 


4.63+ 


0.79 


7.374 


:2.45 




Pope et al. (2008} 


1.910 


0.45 


1.34±0.22 


6.46±1.30 


2.13= 


=0.38 


1.19= 


t0.32 


14.32 


t5.96 


5.90+1.00 




CJoppin et al. (2010} 


2.730 


0.51 


3.04±0.71 


16.28±3.52 


2.55= 


=0.72 










10.98 


±1.88 




harrah et al. (2008 1 


1.690 


0.73 


2.37±0.18 


3.56±0.24 






1.50= 


tO.10 











Brandl et al. (2006} 
O'Dowd et al. (2009 ) 



0.008 
0.092 



(mJy) 

;.52xl0 3 2.10+0.85 
6.08 4.24+0.49 



(10 8 L Q ) 



4.26+1.61 
17.72+1.72 



0.83= 
3.82= 



0.34 
0.36 



1.85= 
5.16= 



0.65 
0.43 



(10 s L e ) 
1.23+0.64 



(10" Le) 
1.18+0.31 



"The values are computed from the stacked spectra. 

"The sample of Sajina et al. (2007} is limited to PAH-rich sources with EW 7 7/nn > 0.8 /jm. 
6 Menedez-D. et al. (2009) designates |Menendez-Delmestre et al. (2009 1. 



Table 7. Average luminosities ratios for different samples. 



Sample 



L 6.2tt 



L l.lii 



L 1.7u 



(io- 3 ) 



vL y (1.4GHz) 

(10 3 ) 



This work 

This work, EW 7 7/lm > 6 yum" 
This work, EW 7 .7m n < 6yum a 
|Sajina et al. (lOOlf 
|Menen dez-Delm estre et al. (2009} 
| Pope et al. (2008| 
Uoppin et al. (2010} 



Earrah et al. (2008 



Brandl et al. (2006 



O'Dowd et al. (2009 



0.28 


1.27 


2.87 


4.53 


1.58 


3.27 


11.70 


11.96 


0.27 


1.22 


3.30 


4.50 


1.36 


3.57 


13.12 


15.46 


0.24 


0.92 


2.78 


3.75 


1.35 


2.16 


8.84 


8.16 


0.20 


1.56 


4.59 


7.70 


1.68 


4.21 


20.84 


7.18 


0.27 


0.78 


2.52 


2.95 


1.17 


1.48 


5.56 


8.86 


0.21 


1.13 


3.03 


5.43 


1.79 


2.27 


10.95 


4.51 


0.19 




6.38 






2.77 


14.83 




0.66 


1.58 




2.37 










0.49 


1.14 


5.13 


2.30 


0.45 


1.77 


3.61 


3.46 


0.24 


0.83 


4.63 


3.44 


0.74 









"The values are computed from the stacked spectra. 

*The sample of Sajina et al. (2007 1 is limited to PAH-rich sources with EW 7 7;lm > 0.8 /jm. 



